* scale_graph_formal.do	JEP 4/10/2018
* Graph formal effects against counterfactual effects
* Formal sector

*************************** Change log ********************************** 

/*
	
*/	

clear 
version 13.1

* Project info
cap project, doinfo
if _rc==198 {
	if c(username)=="J16339" loc master "B:\Col_Minimum_Wage" 
	else if c(username)=="A15858" loc master "B:\Col_Minimum_Wage"
	else if c(os)=="Windows" loc master "C:/Users/jorpp/Dropbox (Brown)/Col_Minimum_Wage"
	else loc master "/home/jperez/AA_Minimum_Wage"
	loc pr=0
}
else {
	local master "`r(pdir)'"
	local doname "`r(dofile)'" 
	* Project calls
	project, relies_on("`master'/Code/Est/graph_rif.do")
	project, uses("`master'/data/source/cfac.dta")
	foreach x in 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 {
		project, uses("`master'/Ster/rif4_new_indcc/reg_emp_bp_tq(1999q1)_bite_salario_mensual_real_50_edom_indcc_p`x'.ster")
	}
	loc pr = 1
}

* Load cfac dataset

use "`master'/data/source/cfac.dta"

drop informality _j

reshape wide wcqe, i(m) j(bicho)

tempfile cq
save `cq'

* Now make formal effects graph and generate vars

cap program drop mex
program define mex
	cap mkdir "`2'/Doc"
	cap mkdir "`2'/Doc/tosubmit2"
	cap mkdir "`2'/Doc/tosubmit2/fig"
	graph export "`2'/Doc/tosubmit2/fig/`1'.pdf", replace
	graph export "`2'/Doc/tosubmit2/fig/`1'.png", replace
	graph export "`2'/Doc/tosubmit2/ILR-Figs/Fig5_`1'.png", replace
	*! epstopdf "`2'/Doc/tosubmit/fig/`1'.eps"
	*erase "`2'/Doc/tosubmit/fig/`1'.eps"
	* project, creates("`master'/Doc/descriptive_figures/density/`1'.pdf")
end	

do "`master'/Code/Est/graph_rif.do"

graph drop _all
glo qtiles "5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90"
set scheme s1color
loc i=1
foreach spec in /*reg_emp_bp reg_emp_bq */ reg_emp_bp {
	if `i'!=1 loc leg "legend(off)"
	if "`spec'"=="reg_emp_bp" loc subt "Bartik price shocks"
	if "`spec'"=="reg_emp_bq" loc subt "Bartik quantity shocks"
	if "`spec'"=="reg_sst_emp_bp" loc subt "Bartik price shocks - City trends"
	if "`spec'"=="reg_emp_bp_ixc" loc subt "Bartik price shocks - City x Industry FE"
	if "`spec'"=="reg_sst_emp_bq" loc subt "Bartik quantity shocks - City trends"
	if "`spec'"=="reg_emp_bq_ixc" loc subt "Bartik quantity shocks - City x Industry FE"
	
	else loc leg "legend(label(2 "UQE") label(1 "95 % CI"))"
	graphrif, master("`master'") ster(Ster) folder(rif4_new_indcc) qtiles($qtiles) spec(`spec') time(tq(1999q1)) rareaopts(fcolor(gs13) lcolor(none) lpattern(dash) ) gen twowayopts( xtitle(Percentile) ytitle(Coefficient) yscale(range(-0.01 0.01)) xscale(range(5 90)) ylabel(-0.01(0.005)0.015, angle(horizontal)) legend(off) name(g`i'), (pci -0.0125 22.66 0.015 23, lcolor(black)) (pci -0.0125 30.98 0.015 31, lcolor(cranberry) lpattern(dash)))   lineopts(lcolor(black) xlab(10(10)90)) 
	loc i=`i'+1
}

gen m=qtile
merge 1:1 m using `cq'

replace wcqe1=0 if wcqe1==.

keep if m>=5

twoway (rarea ul ll qtile, fcolor(gs13) lcolor(none) lpattern(dash)) (line beta qtile, lcolor(black) xlab(10(10)90)) (function y=0, lcolor(gray)  range(5 90)) (line wcqe1 qtile, lpattern(longdash) lcolor(cranberry)) (line wcqe2 qtile, lpattern(dash_dot) lcolor(cranberry)), xtitle(Percentile) ytitle(Coefficient) yscale(range(-0.01 0.025)) xscale(range(5 90)) ylabel(-0.01(0.005)0.025, angle(horizontal)) legend(label(2 "Estimated effect") label(4 "Increase for affected") label(5 "Increase for all below") order(2 4 5)), (pci -0.0125 23 0.025 22.66, lcolor(black)) (pci -0.0125 30.98 0.025 31, lcolor(cranberry) lpattern(dash))

graph play "`master'/Code/Doc/tosubmit2/scale_formal_edit.grec"

mex scale_formal "`master'"

if `pr' project, creates("`master'/Doc/tosubmit2/fig/scale_formal.pdf")


